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Abstract. Arnold, Falk, and Winther [Sm«. Amen Math. Soc. 47 (2010), 281-354] 
recently showed that mixed variational problems, and their numerical approximation 
by mixed methods, could be most completely understood using the ideas and tools of 
Hubert complexes. This led to the development of the Finite Element Exterior Cal- 
culus (FEEC) for a large class of linear elliptic problems. More recently. Hoist and 
Stern [arXiv: 1005.4455, arXiv: 1010.6127] extended the FEEC framework to semi-linear 
problems, and to problems containing variational crimes, allowing for the analysis and 
numerical approximation of linear and nonlinear geometric elliptic partial differential 
equations on Riemannian manifolds of arbitrary spatial dimension, generalizing surface 
finite element approximation theory. In this article, we develop another distinct exten- 
sion to the FEEC, namely to parabolic and hyperbolic evolution systems, allowing for 
the treatment of geometric and other evolution problems. Our approach is to combine the 
recent work on the FEEC for elliptic problems with a classical approach to solving evo- 
lution problems via semi-discrete finite element methods, by viewing solutions to the 
evolution problem as lying in time-parameterized Hilbert spaces (or Bochner spaces). 
Building on classical approaches by Thomee for parabolic problems and Geveci for hy- 
perbolic problems, we establish a priori error estimates for Galerkin FEM approximation 
in the natural Bochner space norms. In particular, we recover the results of Thomee and 
Geveci for two-dimensional domains and the lowest-order mixed method as a special 
case, effectively extending their results to arbitrary spatial dimension and to an entire 
family of mixed methods. We also show how the Hoist and Stern framework allows for 
extensions of these results to certain semi-linear evolution problems. 



Contents 

1 . Introduction 2 

2. Semi-Discrete FEIVI Error Estimates for Evolution Problems 4 

3. The Finite Element Exterior Calculus 6 

4. Abstract Evolution Problems and Bochner Spaces l8 

5. A Pn'on Error Estimates for Parabolic Problems 10 

6. A Priori Error Estimates for Hyperbolic Problems 15 

7. Semi-linear Evolution Problems 19 

8. Concluding Remarks 21 

9. Acknowledgments 21 
Appendix A. Explanation of Semi-Linear Error Estimates 21 
References 23 



Date: February 9, 2012. 

Key words and phrases. FEEC, elliptic equations, evolution equations, nonlinear equations, approxi- 
mation theory, nonlinear approximation, inf-sup conditions, a priori estimates. 

MH was supported in part by NSF Awards 0715146, by DOD/DTRA Award HDTRA-09- 1-0036, and 
by NBCR. 

AG was supported in part by NSF Awai'd 0715146 and by NBCR. 

1 



2 



A. GILLETTE AND M. HOLST 



1. Introduction 

More than two decades of research on linear mixed variational problems, and their 
numerical approximation by mixed methods, recently culminated in the seminal work 
of Arnold, Falk, and Winther Arnold, Falk, and Winther [13J|. They showed that these 
problems could be most completely understood using the ideas and tools of Hilbert com- 
plexes, leading to the development of the Finite Element Exterior Calculus (FEEC) for 
elliptic problems. In two related articles lfT6l[T7l . Hoist and Stem extended the Amold- 
Falk-Winther framework to include variational crimes, allowing for the analysis and 
numerical approximation of linear and nonlinear geometric elliptic partial differential 
equations on Riemannian manifolds of arbitrary spatial dimension, generalizing the ex- 
isting surface finite element approximation theory in several directions. In the current 
article, we extend the FEEC in another direction, namely to parabolic and hyperbolic 
evolution systems. Our approach is to combine the recent work on the FEEC for elliptic 
problems with a classical approach to solving evolution problems using semi-discrete 
finite element methods, by viewing solutions to the evolution problem as lying in time- 
parameterized Hilbert (or Bochner) spaces. Building on classical approaches by Thomee 
for parabolic problems and Geveci for hyperbolic problems, we establish a priori er- 
ror estimates for Galerkin FEM approximation in the natural Bochner space norms. In 
particular, we recover the results of Thomee and Geveci for two-dimensional domains 
and the lowest-order mixed method as a special case, effectively extending their results 
to arbitrary spatial dimension and to an entire family of mixed methods. We also show 
how the Hoist and Stern framework allows for extensions of these results to certain semi- 
linear evolution problems. 

To understand why the finite element exterior calculus (FEEC) has emerged in a natu- 
ral way to become a major mathematical tool in the development of numerical methods 
for PDF, consider the vector Laplacean: 

— Ati = — grad div u + curl curl u, 

and a natural variational formulation: Find u G if (curl; Q) fl i^fo(div; Q) s.t. 

f [(V ■u){V ■v) + {V xu)-{V xv)]dx= [ f-v dx, G H{cm\; n) n i/o(div; fi). 
Jq Jn 

A mixed formulation is a natural alternative: Find (cr, u) e H^{n) x if (curl; fi) s.t. 

[ (ar-u-Vr) dx = 0, \/t E H\n), 
Jq 

[Va ■ V + {V X u) ■ {V X v)] dx = / f-v dx, Wv e H{cm\; n). 
Jci Jq 

Using the standard finite element approach based on the first (non-mixed) formulation 
turns out to not correctly capture either the geometry (in the case of domains with cor- 
ners) or the topology (in the case of non-simply connected domains); there are several de- 
tailed examples in [[3l that explore what goes wrong. A standard finite element approach 
based on the second (mixed) formulation turns out to work extremely well and encoun- 
ters neither of the difficulties of the first approach. Why is this the case? To answer this 
question, consider an arbitrary vector field / G iv^(i^); we have the Helmholtz-Hodge 
orthogonal decomposition into curl-free, divergence-free, and harmonic functions: 



f = Vp + V X q + h. 
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where h is harmonic (divergence- and curl-free). The mixed formulation is essentially 
computing this decomposition for h = and finite element methods are somehow ex- 
ploiting this. 

The relationship can be found by examining the de Rham cohomology of the domain fi. 
The space of harmonic /c-forms is isomorphic to the kth de Rham cohomology group of 
the domain meaning the dimension of the space of harmonic forms indicates the num- 
ber of topological features that can create obstacles to well-posed formulations of elliptic 
problems. A natural question is then: What is an appropriate mathematical framework 
for understanding this abstractly, that will allow for a methodical construction of "good" 
finite element methods for these types of problems? 

The answer turns out to be Hilbert Complexes. Hilbert complexes were originally 
studied in [8| as a way to generalize certain properties of elliptic complexes, particularly 
the Hodge decomposition and other aspects of Hodge theory. A Hilbert complex {W, d) 
consists of a sequence of Hilbert spaces W^, along with closed, densely-defined linear 
maps d'' : <Z — V'^'^^ C W'^^^, possibly unbounded, such that d'' o d''"^ = for 
each k. 

^ yk-l yk ^ ^ . . . 

This Hilbert complex is said to be bounded if d'^ is a bounded linear map from to 
T^k+i ^Qj. k, i.e., (W, d) is a cochain complex in the category of Hilbert spaces. It 
is said to be closed if the image d'^V'^ is closed in for each k. It was shown in [|3l 

that Hilbert complexes are also a convenient abstract setting for mixed variational prob- 
lems and their numerical approximation by mixed finite element methods, providing the 
foundation of a framework called finite element exterior calculus (see also [2]). This line 
of research is the culmination of several decades of work on mixed finite element meth- 
ods and computational electromagnetics OIlOll^IllHl. The most important example of 
a Hilbert complex for our purposes of the FEEC arises from the de Rham complex of 
smooth differential forms on a domain or manifold. 

The main developments in FEEC to date have been for linear (and now semi-linear) 
elliptic problems such as Poisson's equation 

-An = /. 

Our goal here is expand the scope of this analysis to include parabolic linear (and semi- 
linear) equations such as the heat equation, 

{d, - A)u = /, 

and hyperbolic equations such as the wave equation, 

(du - A)u = f. 

The exterior calculus framework treats A as (d + 5)^, where d is the exterior derivative 
operator and 6 its adjoint. The incorporation of the time derivative operation dt into this 
framework, however, has not been previously considered. To remedy this, we develop 
the most natural extension of FEEC theory to evolution problems: a generalization of the 
semi-discrete method often called the 'method of lines.' This approach involves the dis- 
cretization of the spatial part of the differential operator, leaving the time variable contin- 
uous. It can be viewed as introducing a time parameter into the discrete (Hilbert) spaces 
that have been developed for elliptic problems. These parametrized Hilbert spaces are 
naturally treated as Bochner spaces and we will show how Bochner space theory makes 
obvious the well-posedness of the problems we consider. Moreover, the accompanying 
Bochner space norms, when coupled with FEEC notation for Hilbert complexes, provide 
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a clear and consistent notation for bounding errors in mixed methods accumulated over 
a finite time interval. 

We note that there is another approach to solving evolution problems with finite ele- 
ments, namely using a complete discretization of space-time. This tactic allows for the 
dynamical change of the underlying discrete approximation spaces in both space and 
time. Such an approach gives rise to space-time adaptivity, and is potentially the most 
flexible and powerful approach to the numerical treatment of parabolic and hyperbolic 
evolution problems. This approach, which we will consider in a second article, is most 
naturally formulated using geometric calculus, a well-studied mathematical structure for 
time-dependent problems. In the current article, we focus on extending FEEC to semi- 
discrete methods using Bochner norm estimates for the method-of-lines approach. 

Finally, we note that the work presented here was developed simultaneously and inde- 
pendently from a related project by Arnold and Chen [JJ. While our focus is extending 
the scalar Hodge-Laplacian to both linear and semi-linear evolution problems, with or 
without geometric variational crimes, their focus is on extending FEEC theory to gener- 
alized Hodge-Laplacian parabolic problems. The pairing of these two results will lead to 
further insight on the subject in a variety of directions. 

Outline of the paper. The remainder of the paper is structured as follows. In Sec- 
tion [21 we review the classical semi-discrete mixed finite element method error estimates 
for parabolic problems due to Thomee [|27] | and others for hyperbolic problems due to 
Geveci [12J. In Section [3l we give a very brief overview the Finite Element Exterior Cal- 
culus and recall some relevant results. In Section |4l we formulate abstract parabolic and 
hyperbolic problems in Bochner spaces and state some standard results on the existence 
and uniqueness of strong and weak solutions. In Section [51 we combine the classical 
approach to semi-discrete methods with modern FEEC theory to establish some basic 
a priori error estimates for Galerkin mixed finite element methods for parabolic prob- 
lems. The main result is Theorem 15.21 which exploits the FEEC framework to obtain a 
classification of spatial finite element spaces that give optimal order convergence rates in 
Bochner norms. In Section [6l we carry out a similar analysis for hyperbolic problems, 
resulting in the error estimate given in Theorem 16. 2[ Our results recover the estimates of 
Thomee and Geveci for two-dimensional domains and the lowest-order mixed method 
as a special case, effectively extending their results to arbitrary spatial dimension and to 
an entire family of mixed methods. In Section [71 we employ the results of Hoist and 
Stern to extend our parabolic estimates to a class of semi-linear evolution PDF. Finally, 
in Section m we draw some conclusions and make some remarks on future directions. 

2. Semi-Discrete FEM Error Estimates for Evolution Problems 

We begin by reviewing semi-discrete finite element methods and their a priori error 
estimates for parabolic and hyperbolic PDE systems. We focus in each case on a rela- 
tively simple, well-studied system of interest to modeling communities, namely, the heat 
equation (parabolic) and the wave equation (hyperbolic). The heat equation is: find 
u{x, t) such that 



We review the approach to Galerkin methods for this problem as presented in Thomee [[271 
for domains VL C M^. His approach is based on an earlier work with Johnson [18] and 
builds upon prior analysis of elliptic projection [7J. Let a = Vu and define the mixed. 



Ut — Au = f inQ, for t > 

u = on d^l, for t > with u{ 



0) = ^ in n. 



(2.1) 
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weak form problem: Find {u, a) E L'^ x H{div) such that 

K»-(diV0r,(/)) =(/,0), t>0 

(a, + (m, divw) =0, VwGi/(div), t > 0, u{0) = g. ^ ' 

The semi-discrete problem is then to find [uh, cfh) E Sh x Hh C L"^ x H{div) such that 

{uh,tAh) - {divah,(l)h) ={f,(f)h), ^(t)h^Sh, t>0 
{(7h,uJh) + {uh,diwuh) =0, \/uhEHh, t>0, Uh{0) = gn- 

where gh is an approximation of g in Sh- With bases for Sh and Hh, the matrix form of 
the discrete problem is 

AUt - bj: = F, 

B^U + DS = 0, for t > 0, U{0) given, 

where U and S are vectors corresponding to Uh and ah- It is easily seen that the matrices 
A and D are positive definite. Eliminating S, we have the system of ODEs 

AUt + BD-^B^U = F, for t > 0, U{0) given, 

which by standard results in ODE theory has a unique solution. 

Thomee uses discontinuous linear elements for Sh and piecewise quadratic elements 
for Hh- He defines the solution operator Th : L?' Sh given by T/j/ = Uh for the 
corresponding elliptic problem. Then gh can be defined by 

gh := Rhg ■= -ThAg. 
For gh = Rhg and t > 0, Thomee derives the estimates 

\\uh{t) - u{t) 1 < ch^ (^\\u{t) + \Wt\\md^, (2.4) 



Note that these estimates are for a fixed time value t and restricted to a particular choice 
of finite elements in 2D. 

We now turn to the wave equation: find u{x, t) such that 

Utt — An = / mVt, for t > 0, 

u = Q on dVL, for t > with 0) = uq in (2.6) 

and 0) = til in f2 

There are two approaches to defining a mixed weak form of this problem. The first is 
very similar to the parabolic case: given /, uq, and ui, find (u, a) such that 

K,0)-(diva,(/)) =(/,0), t>0, 

(a, + (u, div w) =0, Vw G //(div), t > 0, 

u(0) = Uq. 

ut{0) = ui. 



(2.7) 



It is difficult to derive estimates for the numerical approximation of (12.71) akin to those 
found in the parabolic case due to the second derivatives appearing in the formulation. 
Some estimates attempts along these lines for C have been given by Baker [4] and 
Cowsar, Dupont and Wheeler ifTOlfTTIl . 
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For the purpose of extending the FEEC framework, we find the 'velocity-stress' for- 
mulation of the problem and the resuhs of Geveci lfT2ll to be more useful. This formula- 
tion solves for fi := ut instead of u: given /, uq, and ui, find cr) G x i/(div) such 
that 

(/ii»-(diva,0) =(/,0), t>0, 
(crt,a;) + (/i,div w) =0, VwGif(div), t > 0, 

(2-8) 

^(0) = Ml, 
a(0) = Vmq. 

The semi-discrete problem is then to find (/i^, ah) E Sh x Hh C x iJ(div) such that 

iHh,t, (l>h) - (diy CTh, (l)h) =(/,0h), \/(l)hESh, t>0, 
{ah,t, i^h) + (/"h, div Uh) =0, \/ Uh e Hh, t > 0, 

(2,9) 

/x/j(0) = ni,,,, 
a,,(0) = {Vuo)h, 

where ui^h is an approximation of ui in and (Vuo)h is an approximation of Vuq. 
Again, bases for Sh and Hh reduce the discrete problem to a matrix formulation: 

AWt -BE = F, 

B'^W + DY.t = 0, for t > 0, W{0), S(0) given, 

where W and S are vectors corresponding to fih and ah and A, D are symmetric, positive 
definite matrices. As Geveci [fT2L p. 248] explains, this can be reduced to a single iterative 
system of the form 

{D + eB^A-'B)J:"^' = G, 

where k denotes the time step in an implicit Euler time-differencing scheme. 

To derive an error estimate for the velocity-stress discretization, Geveci states the need 
for projection operators from i7(div) to Hh and from to Sh satisfying certain approx- 
imation properties. He explains that such operators exists for a variety of finite element 
spaces in M^, e.g. the Raviart-Thomas spaces f,22,l . allowing the following result. For 
1 < s < r with r > 2, 

\\fih{t) - Kt)\\L2 + \Wh{t) - < c(||mi - Ui,h\\ + IIVmo - {Vuo)h\\) + 

+ch^ (^\\u,\l + \\Vuo\l + l\\\f^t{r)\l + \\atir)\ldr 

(2.10) 

Like estimates (12.41 ) and (12.51) for the parabolic problem, (12.101 ) says that the approxi- 
mation error can be controlled in norm at any time t by the H'^ norm of the initial 
conditions plus the accumulated norm of the variables up to time t. It is these types 
of estimates that the FEEC framework can refine, simplify, and generalize to arbitrary 
spatial dimension n. 

3. The Finite Element Exterior Calculus 

The finite element exterior calculus (FEEC) provides an elegant mathematical frame- 
work for deriving error estimates for a large class of elliptic PDF. We now give a brief 
overview of the notation and certain main results from FEEC which are relevant to this 
paper; we refer the reader to the seminal papers of Arnold, Falk, and Winther ^ [3]| for 
additional explanation. 
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Let ^7 be a bounded n-manifold embedded in M" and assume has a piecewise 
smooth, Lipschitz boundary. The space of -bounded continuous differential k-forms 
on i7 is given by 

L'^K^{n) ■= |^a/(ix/ G A\n) : ai G L\Q) V/| , 

where / ranges over all strictly increasing sequences of k indices chosen from {!,...,«}. 
The exterior derivative operator dk : A''{n) — )• A''+^(r2) acts on these spaces to form a 
Hilbert complex (L'^A,d). The associated domain complex is the sequence of spaces 
HA'' := domain(dfc) C L^A'^(i7), commonly called the deRham complex: 

. //AO ^ i/Ai ^ . . . ^ //A- . 0. 

The norm on each space is the graph norm associated to d, i.e. 

{u,v)Hk(^n) ■■= {u,v)Ak(n) + (4m, 4v)Afc+i(f7)- 

We note that in any dimension n, the beginning and end of the deRham complex can 
be understood in terms of traditional Sobolev spaces and differential operators: 

. H\Q) ^ . . . ^ L\Q) 

A major conclusion of FEEC is that stable finite element methods for elliptic PDE must 
seek solutions in finite dimensional subspaces A| C HA'' that satisfy certain key ap- 
proximation properties. First, the subspaces should form a subcomplex of the deR- 
ham complex, meaning dA^ C A^^^. Second, A^ should have sufficient approximation 
that upper bounds on inf„gy\^fc \ \u — v\\jj^k can be ensured for some or all u G HA''. 
Third, there must exist bounded cochain projections vr^ : HA'' A^ which are in- 
variant on A^, commute with the exterior derivative operators, and provide a bound 

In the context of the deRham complex, all these properties are shown to be provided 
for by two canonical classes of piecewise degree r polynomials associated to a simplicial 
mesh T of Let Vr denote polynomials in n variables of degree at most r and "H,. C Vr 
the subspace of homogeneous polynomials. The first class, denoted VrA''{T), consists 
of all /c-forms with coefficients belonging to Vr on each n-simplex of T. The second 
class, denoted 'P^A^(T), interleaves with the first class, i.e. 

Vr-iA\T) C V;A\T) C VrA\T). 

To define V~A''{T), first define X be the vector field on such that X{x) is the vector 
based at x G M" that points opposite to the origin with length Define A''{T) := 
VrA'' © K'Hr-iA''^^ , a direct sum, where k is defined by contraction with X. The map k 
is called the Koszul differential and gives rise to the Koszul complex. This is elaborated 
upon in detail in the work of Arnold, Falk and Winther e.g. [3, p. 328]. 

For n = 3, we have the following correspondences between the FEEC notation of 
finite element spaces and traditional element spaces. 

Vr+iA^{T) = Nedelec 2nd-kind (div) elements of degree < r + 1 (see (HI) 

Vr+iA'^{T) = Nedelec Ist-kind if (div) elements of order r (see [|20ll ) 

Vr+iA^^{T) = discontinuous elements of degree < r 

VrA?{T) = discontinuous elements of degree < r 
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Hence, in the case of the deRham complex, FEEC recovers the well-known finite ele- 
ment spaces while at the same time describing their generalization to arbitrary spatial 
dimensions. 

The last piece of FEEC used in this work is the existence of smoothed projection 
operators 

Til : L^A'^ ^ Kl where A^, e {n.A'^(r), V;X\r)]. (3.1) 
The following theorem asserts the key properties of these operators. 

Theorem 3.1 (S Theorem 5.9). 

(1) Let Kl be one of the spaces V;^^k''{T) or, if r > 1, Vrh-^{T). Then tt^ is a 
projection onto A^ and satisfies 

for < s < r + 1. Moreover, for all uj G LP'A^iyi), nf^u — )■ u in L"^ as h ^ 0. 

(2) Let A-l be one of the spaces VrA^{T) or V; A\T) with r > 1. Then 

for < s < r. 

(3) LetKl-^ e {Vr+iA^~\r),V;^^K^-\T)] and = V;^^K\T) or, ifr > 0, 
VrA\T). Then dvr^^^ = n^d. 

An explicit construction of these operators can be found in the papers of Arnold, Falk, 
and Winther [Em. 

4. Abstract Evolution Problems and Bochner Spaces 

We now cast parabolic and hyperbolic problems into the abstract framework of Bochner 
spaces, i.e. spaces of functions whose values lie in a Banach space. This approach to an- 
alyzing parabolic and hyperbolic evolution PDE can be found for example in [|26l page 
66] and [23 J. 

Let X be a Banach space and / := (0, T) an interval with closure / := [0,T]. Define 

C{i, X) := {u : i ^ X I u bounded and continuous}. 

Equip this space with the norm 

ll^llc(i,x) := sup||tx(t)||^. 
tei 

The Bochner space L^{I, X) is then defined to be the completion of C{I, X) with re- 
spect to the norm 

, i/p 

MLPii,x) '■= ( / \Ht)\\xdt 



The space H^{I, X) has an analogous norm 




|2 



u\\m{i,x) ■= I L\Ht)\\x + 



d . . 




We now consider abstractions of the two main types of evolution PDE. 

Let H and V be real, separable Banach spaces such that V is continuously and densely 
embedded in H. This provides a Gelfand triple, also called a 'rigged Hilbert space' : 

V CH CV*, 



FEEC FOR EVOLUTION PROBLEMS 9 

where H is also continuously and densely embedded in V*. It should be noted that the 
isomorphism between V and V* is in general not the same as the composition of the 
inclusion mappings. 

Let (-, ■) denote the inner product on H as well as the natural pairing between V* and 
V. Given A{t) e C{V, V*) depending continuously on t e /, define a quadratic form 

a{t,u,v) := -{A{t)u,v), (4.1) 

for {t,u,v) G W X V X V. Assume that a satisfies the coercivity condition 

a{t,u,u) > Ci \ \u\\y — Ci IImII^ , (4.2) 

with ci, C2 constants independent oftEl. Consider the abstract parabolic problem 

ut = Ait)u + fit), (4.3) 
^(0) = Uo, (4.4) 

and the abstract hyperbolic problem 

u = A{t)u + f{t), (4.5) 

u{0) = Uo, (4.6) 

u{0) = Ml. (4.7) 

These abstract formulations are well-posed in the following sense. 

Theorem 4.1 (Existence of Unique Solution to the Abstract Parabolic Problem). Let 

f e L'^{I, V*) and uq e H, let a{-, ■, ■) be as in (PTTI) . and let (|431) hold. Then the 
abstract parabolic problem ( 14.31) has a unique solution 

u e L'^{i,v)nH\i, V*). 

Moreover, the Sobolev embedding theorem implies u E C{I, V*), allowing an interpre- 
tation of the initial condition u{0) = uo. 

Proof See [lH, page 382. □ 

The analogous result for the abstract hyperbolic case will require two additional con- 
ditions: 

a{t,u,v) = a{t,v,u), 'i u,v eV, (4.8) 
Ae C^(!,C{y,V*)). (4.9) 

We now have 

Theorem 4.2 (Existence of Unique Weak Solution to the Abstract Hyperbolic Problem). 
Given f G L^(J, H), uq G V, and ui G H. Let a(-, -, ■) be as in (|4l1) . and let (|42|) . (l48l) . 
and (14.91 ) hold. Then the abstract hyperbolic problem ( 14.51) has a unique weak solution 

u G C(7, y) nc^(7, H). 

Proof See BH, page 389. □ 

These results give a short and elegant way to prove that a wide class of PDE problems 
amenable to finite element methods are well-posed. 



10 



A. GILLETTE AND M. HOLST 



5. A Priori ERROR ESTIMATES FOR PARABOLIC PROBLEMS 

We extend Thomee's error estimates from Section [21 to the broader class of elements 
and arbitrary spatial dimension allowed by FEEC using the abstract framework estab- 
lished in Section m Fix an interval / = [0,T] and define the Bochner differential 
A;-form space by parametrizing the deRham complex and its domain complex. For 
these spaces, we use the notation 

X'' := L\I,L^ A\n)) and HX'' := L\I , HA^Q)). (5.1) 

We first observe that the parabolic problem (12.11) is a special case of the abstract parabolic 
problem (|4.3I) - (|4.4I) . Thus by Theorem 14. 1[ it has a unique solution on (0, T) with a C°- 
continuous extension to [0, T]. We may therefore denote the interval of solution to be 
/ := [0, T] without confusion. For i7 c M", define the Bochner mixed weak parabolic 
problem: Given / and g, find (m, a) G i/X" x HX^'^^ such that 

(wt,0)-(diva,0) =(/,0), V(/)G//A", tel, 

{a, u) + {u, div oj) =0, \/ u e //A"~\ t e I, (5.2) 

u{0) =g. 

The semi-discrete Bochner parabolic problem is: Find {ut, at) e x A^~^ such that 
{uh,t, (ph) - {diy (Th, (f)h) ={f,(ph), ^(t)h^K^ tel, 
{ah, Uh) + {uh, div Uh) =0, ^ Uh e A^-^ t e /, (5.3) 

Uhio) = gh- 

where E A^ is an approximation of g. Define gh to be the solution to the elliptic 
problem with load data — A^', i.e. 

{div ah Ah) + {Ag,(t)h) =0 V0;,gAJJ, 

(5-4) 

{ah, Uh) + {gh, div Uh) =0, WuJhE Al . 

Viewed another way, gh is the image of the solution operator Th : L"^ ^ A^ for the elliptic 
problem, i.e. gh = —ThAg. It can be shown, in the same manner presented in ETl . that 
Th is positive definite on A^ implying that the semi-discrete Bochner problem has a 
unique solution. 

Elliptic projection, an idea dating back to Wheeler [|29ll , can be carried out for any 
fixed time value as we now discuss. For any Iq e I, define the time-ignorant discrete 
elliptic problem: find {uh, G A^ x AJJ~^ such that 

{div ah, 4>h) + {-Au{to),(f)h) =0, \f(/)hEAl, 

{dh, Uh) + {uh, div Uh) =0, WuhE Al'^, (5.5) 

Uh{^) = gh- 

Note that the u appearing in the first equation of (15.51) is the solution to the continuous 
problem (15.21) . Thus, we can view ah and Uh as functions of t with the understanding 
that they are defined for each t value by (15.51) alone; no continuity with respect to t is 
required, hence the moniker 'time-ignorant.' 

For ease of notation, and in keeping with Thomee, define the error functions 

p{t) := Uh{t) - u{t), 

e{t) ■.= Uh{t)~Uh{t), 

e{t) := ah{t) - ah{t). 
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We now prove a lemma which will aid in our subsequent analysis; the result for the linear 
parabolic problem appears as part of the proof of Thomee Il27l Theorem 17.2]. 



Lemma 5.1 (Thomee 11271 '). The error functions satisfy the semi-discrete formulation 

(0^,0/,) - (dive, 0/,) = -(pi, V0/,gA;J, tel, 



(e, uoh) + {e, div Uh) =0, M Uh e Kl ^ , t e 1. 



(5.6) 

h 1 ^ ' 



Proof. The second equation is immediate from the second equations in (15.31) and (|5.5I) . 
The first equation can be written out as 

(t^h) - (div ah, (j)h) + (div ah, (ph) - {uh,u <Ph) = {ut, 0h) - (uh^t, <Ph) 
which is reduced as follows: 

{uh,t, (t)h) - (div ah, (ph) + (div ah, (ph) = {ut, <ph) cancel like terms 

iuh,u <Ph) - (div ah, (ph) = -(An, (ph) + (ut, <ph) by (1531) 

if,(ph) = -i^u,(Ph) + iut,<ph) by(l53]) 

This says that the continuous problem ut — Au = f should hold in a weak sense when 
tested against any of the functions in A^. This is guaranteed to be true since we chose 
C A" = L^. Thus, the error equations hold as stated. □ 

The following theorem says that if A^ and A^^^ are chosen according to the FEEC 
framework, then error estimates akin to Thomee's can be obtained. 

Theorem 5.2. Fix I := [0,T]. Suppose {u,a) is the solution to ( 15.21) such that the 
regularity estimate Q 

+ + IkWIIn^+i + WdcrmHs < c||/(t)||^. (5.7) 

holds for < s < Smax <^nd t E I. Choose finite element spaces 

K-' = { or }, Ki = v;^,K-{r) (=nA"(r)) 
n+iA"-Hr) 

Then for < s < Smax, Qh defined by ( 15.41) . and {uh,ah) the solution to ( 15.31) . the 
following error estimates hold: 



ch^^' (\\f\yii,Hs) + ^m\L^ii,H^)) forr >0,ifs<r-l 

(5.8) 



^The regularity assumption ( 15.71 ) comes from Arnold, Falk, and Winther 13] p. 342]. Note that du = 
since m e A"; we leave the associated term only to show consistency with ||3l. 



12 



A. GILLETTE AND M. HOLST 

+ V^II/*II 

ifr = a s = a a;^'^ = A"-i(r) 

ifr = 0, s < 1, = PiA"-i(r) 



C /l 



L2(/ 



for r>0, ifs<r — 1 



(5.9) 



L2{I,H^) + ^ ll/tllL2(/_L2) 



|div((Tft, - cr)||i2^n < < 



i/r = 0, s < 1 

for r>0, ifs<r — 1 



(5.10) 



Proof. To simplify notation, we will often use 1 1 • 1 1 to mean 1 1 ■ | |/^2(j^) . We adapt the proof 
technique of the corresponding theorem from Thomee ||27l Theorem 17.2] to our setting. 

Observe that (15.51) is exactly the k = n case of the Hodge-Laplacian problem analyzed 
by Arnold, Falk and Winther [3J and the hypotheses here match their hypotheses. We 
can thus use a triangle inequality argument for each estimate, e.g. 

\\uit) - uhimL^ < \Hi) - MmL-^ + Wuhit) - Mmi^ = iipwil^ + ii^wiIl^ 

(5.11) 

The first term will be bounded using the estimates from Q and the second by the tech- 
niques from Thomee l|27l . We can bound \ \p{t) \ |^2 easily by applying the Arnold, Falk, 
Winther estimate [3, p. 342], which gives 



IIpWII 



L2 



< 



chwmw^. 



if r = 

if s < r — 1, for r > 



(5.12) 



Bounding ||6'(t)||^2 is more subtle. Set (ph '■= and Uh := e in (15.61) . Adding the 
equations yields 

~\\0\\' + \\e\\' = -{Pt,0), tel 

We use a technique from Thomee ETl p. 8] to derive an estimate for \ \9(t) \ \ . Since 
may not be differentiable when 9 = 0, introduce a constant 5 > and observe that 



the last step following by Cauchy-Schwarz. Since 1 16'| | < 



-- e r < 



\pt\ 



1 d 

^ + 5^)^/^ we have that 



\Pt\ 



Note that 9(0) = Uh{0) - Uh{0) =9h-9h = 0. Thus 

d 
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Using the bounds on \ \p{t) \ \ from (I5.12L we get 

ch fwfmk^M ifr = 0, 
< < (5.13) 



c/i^^' / \\ft{^)\\H'M forr > 0,if s < r- 1. 



We can now assemble estimate (15.81) by collecting our results. We show the technique of 
the case r = as the other case employs identical analysis. 



Uh-u\\^2xr.= [ / \\uh{t)-u{t)\\ dt 



1/2 

2 



< ( / mm + mmfdt 



1/2 



<ch \ I (ii/(t)ii+ / \m)\\L^Mj dt"^ 
<ch \ I 2|ii/(t)ii^ + ( / \m)\ydi] ]dt 

Roll the 2 into the constant c and observe that the inner integral is maximal when t = T. 
Thus, 

<ch( I \\fm' + \\ft\\hiiM)dt 



1/2 

1/2 



= c/i(ll/lli^(/,L^)+^ll/*llW, 
= c/i(||/||^.(,,i2)+V^||/i|L.(,,^.)). 



We now turn to (15.91 ). i.e. an error bound for the approximation of a. We use the 
same technique of bounding \ \a{t) — ah{t) \ \ by the corresponding Arnold, Falk, Winther 
estimate and ||5"/i(t) — cr/i(t)|| (= ||£(^)||) by a modification of (15.61 ). First, observe that 
the Arnold, Falk, Winther estimate (31 p. 342] gives 



|a(t)-a,(t)||<c/i^+^||/(t)||^.,if 



s<r + l, = P,+iA"-i(r) 



5<r, ^l-' = v;^,^-~\T) 

(5.14) 

To bound ||e(t)||, differentiate the second equation of (15.61 ) with respect to t and set 
(j)h := 2et, ujh := 2£ yielding 

(^i,2^i)-(div£,2^^i) =-(pi,2^0, 
{et,2e) + {9t,div e) =0. 
Adding the equations and converting to norms, we have the bound 

^||e||' + 2||^i||' = -2(pi,^,)<IIPtll' + ll^t||', 
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by Cauchy-Schwarz and the AM-GM inequality. Note that since 6{0) = 0, we have 
£(0) = by the second equation of (I5.6I ). Thus 

\m\f = £ ^ Ms)\fds < ^^Ms)\f + mifds < l^iptifds. 

As before, we use (15.121) to derive 



f / rt \ V2 

Chi / mmUi] ifr = 0, 







1/2 (5.15) 



m\\ < 



We assemble estimate (15.91 ) by collecting these results in the same fashion we did for 
(|5.8I) . In the r = cases, we have 



^^^^'Uo "^*^^''"^°^^) forr > 0,if s < r- 1 



kft - a||^23en-i = ( y \\ah{t) - a{t)\\ dt\ 



1/2 

< ( / {Mt)-a,m + \Ht)\\fdi 




2 \ 

<ch{ I 2(h'^\\m\\i. + l\\f,mUiyt^ ' . 

Rolling the 2 into the constant c and again noting that the inner integral is maximal when 
t = T, we recover the first two estimates of (15.91) : 

/ .T X 1/2 

\Wh - < ch ^ h'^wfmi.^ + dtj 

When r > 0, (15.141) requires s<rors<r+lto obtain optimal convergence rates on 
the first term of the right side while (15.151) requires s < r — 1 to obtain optimal rates on 
the second term. Thus the hypothesis s < r — 1 implies both (15.141 ) and (15.151 ); the last 



estimate of (15.91 ) then follows by identical analysis to the first two cases. 

Finally, we turn to estimate (|5.10l) and follow the same technique. Since div is the 
exterior operator d between the spaces A"^^ and A", we have from Arnold, Falk, and 
Winther [3,p. 342] that 

||div(a(t) - < ch^Wfrnns , if s < r + 1. (5.16) 

To bound ||div £|| set Wh ■= £ in (|5.6I) and take the derivative with respect to t. This 
yields 
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Note that div 5 G since the discrete spaces are chosen to satisfy the relationship 
di\A^~^ C AJJ. Thus, we can set (ph ■= div e in in (15.61) and substitute to get 

llell^ + ||div e\f = (pt,div e). 

at 

By Cauchy-Schwarz, we have | |div £| |^ < \\pt\ \ | |div e| | and hence 



|div e\ \ < \\pt\ \ ■ 



Again, we use (15.121) to get 



ch\\ft\\r2 ifr = 
|dive(t)ll<<! ^ , , . (5.17) 

ch llftW^s if s < r- l,forr > 

The estimate (15.101) follows by combining this with (|5.161) . □ 

6. A Priori ERROR ESTIMATES FOR HYPERBOLIC PROBLEMS 

We now analyze hyperbolic problems using the same Bochner differential form 
spaces L^X*' introduced in (15.11) . Observe that the wave equation (12.61) is a special case 
of the abstract hyperbolic problem (I4.5I) - (I4.7I) . Thus by Theorem 14. 2[ it has a unique 
solution on the interval / = [0,T]. For Q c M", define the Bochner velocity-stress 
mixed weak formulation: given /, uq, and ui, find (/i, a) G HX"^ x ifX"~^ such that 

(/ii,0)-(diva,0) =(/,0), W(t>eHA\ tel, 

{at, uj) + {p, div u) =0, Va;Gi/A"-\ tel, 

(n\ (6-1) 

/i(0) = Ml, 

a(0) =Vuo, 

where p = ut as. in (12.81) . The semi-discrete Bochner hyperbolic problem is: Find 

{ph, (^h) e A;J X A^-^ such that 



0/.) - (div ah, 4>h) = (/, (ph), ^ (ph^Al, te I, 

{ah,t, uJh) + (Ph, div Uh) =0, V a;^ G A;^~\ t G /, 

/ift,(0) = Ml,,,, 

a(0) =(Vmo)/». 



(6.2) 



The following lemma recasts notions about properties of projection operators and dis- 
crete spaces assumed by Geveci in [[T2l into the language of FEEC. 

Lemma 6.1. Choose finite element spaces 

K~' = { or ), Ai = v;^,A^{r) (=n.A"(r)). 




Then the projection operators from di.il) satisfy 

{div Txl-^uj, (ph) = {div co,(f)h), V w G H{div), (ph G A^, (6.3) 
«0, div Uh) = (0, div Uh), \/ (j) e , Uh e Al'^ , (6.4) 
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with the approximation properties 

\\n]l-^u-u\\ ^2^„_i <ch'\\u\\ j^s^^-i , 

forO< s<r + 2, ifA"?-^ = Vr+iA''^\r), or 

(6.5) 

forO<s<r + l, ifA^-' = V;^^A^-\T), 
<c/i^||0||^.A„, /orO<s<r + l. (6.6) 

Proof. We denote div by d since it is the exterior derivative operator in this context. By 
Theorem l3.1l (l). vr^^^ and ttJ^ are projections onto Hilbert spaces and hence self-adjoint. 
By Theorem [3TT] (3), we have that 

{dTTl~'^w,(f)h) = {n^du,(f)h) = {du^Tilcph) = {doj,(f)h), 

establishing (|6.3I) . A similar trick establishes (16.41) . Estimate (|6.5I) follows directly from 
Theorem l3.1l (l). (Note that Theorem 13. 11 (1) is stated for the case P^A""^ while here we 
have thereby allowing for the higher bound on s in this case.) Finally, since 

tt'^H — fi = dt{T^^u — u). Theorem [3711 (1) also implies (16.61 ). □ 

Theorem 6.2. Fix I := [0, T]. Choose finite element spaces 

AV={ or I, A^ = p-^,A"(r) (=n.A"(r)). 

Then for (/i/j, ah) the solution to A6.2\l . the following error estimate holds: 

(yfE^ + h'{VTE2 + E3 + VtE^)^ , (6.7) 

where 

Ei = \\ 

Ui — U\^y}\jj2. + ll(Vito) — (V'Uo)/i| |£^2 {error due to discretization of initial data) 
-E'2 = 1 1 "ifi 1 1 + 1 1 Vuo 1 1 s ( regularity of initial data ) 

-^3 = 1 1""* I H=) + 1 1^1 jys) {regularity of continuous solution to d<5.i|) ) 

^4= II ^**llL2(7^j:fs) + 1 1*^* I Il2(j ji^s) {dt of E^ accumulated over I). 

Proof. We adapt the proof technique from Geveci [fT2ll to the FEEC setting. Define \E' := 
p^n ^ y^n-i \vith finite dimensional subspace '■= x A)J^^. Denote the components 
of an element G ^ by {0^, Wj}. The inner product and norm on ^ are 

(6,6)* := (0i,02)l2 + (wi,a;2)L2 and := ^(60*- 

Define a skew- symmetric bilinear form a : \E' x — )■ R by 

a(6,6) := -(div Ui,(I)2)l2 + (0i, div a;2)L2. 

Let tjj := {fi, a) G \E' be the solution to (16.11) and let ^ := (0, a;) G \E' be arbitrary. Then 
adding the equations of (16.11) yields 

(^t,O* + a(V^,O = (/,0)L2 V^Gv]/. (6.8) 
Similarly, from (16.21) we get 

{A,t, Chh + aii^h, Ch) = (/, h)L^ V a e (6.9) 

Define a projection operator tt/i : ^ — )• using the bounded cochain projections from 
(l3ni) via 

7r,e:=«0,<-ia;). 
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By Lemma [6m we have that 

Further, since tt/j only affects the spatial variables, it commutes with the time derivative 
operator, i.e. 

Combining these identities with (16.81) . we have that 

{dt'Kh^{t),ihh + a{T:h^{t),ih) = {f{t),(f)h)L^ + {'iThdMt)-dt^{t),^hh, Va e 

(6.10) 

Now define the error function 

Shit) ■■= TThlpit) - Tphit). 

Subtracting (|6^ from (16.101) yields 

{dteh{t),^hh + a{£h{t),Ch) = {'r^hdMt)-dt(j){t),Chh, V^e^h. (6.11) 
Define an operator Eh : — ^ ^/iby 

Letting Ph : ^ "^h denote projection with respect to (-, ■)^, we can re-write (16.1 II) as 

dtShit) + EhShit) = Ph {-Khdtm - dt(t){t)) (6.12) 
Now, since a is skew-symmetric, Eh is skew-adjoint, i.e. 

Hence, —Eh will generate the one-parameter unitary group e~*^'' . Therefore we can turn 
(16.121) into 

eh{t) = e-'-^eh{Q)+ [ e-'^'-^^~'^Ph{'KhdA{r) ~ dr(t){r)) dr. (6.13) 
Jo 

Since e~*^'' is unitary, it preserves ^-norm, i.e. ||e~*^'"0||^ = HV'IU if) E . 

Thus (16.121) becomes 

\\eh{t)%<\\^hmU+ [ \\7^hdr(pir)-dr(p{r)\\^dT. (6.14) 

Jo 
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Using (16.51 ) and (16.61) from Lemma [6n we start building up the main result. We use the 
notation / < to mean f < eg where c is some constant independent of h and T. 

\2 



ll^h ~" fJ'\\l2j^n + \ \(^h ~ ^IIl2^ 

<2 Wfih- fi\\l2 + \\(Jh- aWlidt 
Jo 

^ / {\\f^h-Ki^\\L2 + \\'^h-K''^ 

Jo 



2 I I —1 I 1 2 



i-T 

/ \\^h{t)\\l2dt+ ||vr>-^||i2i-„ + ||<"V-or||^2^„_iC?i 
Jo 



< / \\eH{mi^ + h'n\m\l^ + h'nwmisdt 

Jo 

< [ Wshiml.dt + h^^El (6.15) 
Jo 

To bound | |£/t(^) 1 1^2, we use (16.141 ) and the same results from Lemma [6711 to get 

W^him < \Kui-ui,h\\ + ||<-'(V?io) - {Vuo)h\\ 

+ [ \KdrKr)-dMT)\\ + \\7r]l-'dMr)-dMr)\\dT 
Jo 

\ui -ui^hW + II(Vmo) - i^uo)h\\ + /i' + IIVmoII^O 
+ (^l\\dMr)\\Hs + \\draiT)\\^. dry 



< 



Integrating this from to T yields 

\\e,m'dt<TE',+h'^TEi 



+ h^' (^j\\dMr)%s + \\dMr)\\HsdT^ dt 

< TE, + K'^TEl + h^^ j^^ (p./^IIl2(,,^.) + lia.^11.2,, „s^^dt 



<TEi + h^'TEl + h^'TEl (6.16) 
Combining (16.151) and (16.161) yields 

{\W - yu||L23-n + \\(Th- (tWl^x^-.? < c [TEI + h^'TEl + h^' EI + /i2^T£;2) 

< c (Vt^i + /i^ (^Vr^a + ^3 + VtE^ ) ' 

Taking the square root of both sides completes the proof. □ 

Theorem 16 .21 implies that the error in ut (i.e. ji) and the error in a due to discretization 
of the continuous problem is bounded a priori over [0, T] by terms depending only on 
the regularity of the initial data and the continuous solution. Since the interval is finite, 
this implies control over the error in u as well. These results imply that the finite element 
choices allowed by FEEC theory for elliptic problems also result in stable semi-discrete 
finite element methods for hyperbolic problems. 
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7. Semi-linear Evolution Problems 

We now show how the techniques developed above can be extended to certain types 
of non-linear evolution problems. Consider the semi-linear heat equation: Find u{x, t) 
such that 

-ut - A-u + F{u) = / in for t > 

u = Q ondVt, fort>0 with 0) = ^ in ^ ' 

where F is some non-linear operator on L'^iVt). The existence and uniqueness of solu- 
tions to instances of this problem have been studied extensively [[TSl [T3l l24l[30l as have 
finite element methods for the approximation of its solution ll28l l9l l25ll27l[T9l . 
We focus here on the case where F satisfies a Lipschitz condition 

\\F{v)- F{w)\\^,<C\\v-w\\^2, \fv,weL\Q). (7.2) 

This condition, or a weaker locally Lipschitz condition 10, is assumed by Hoist and 
Stem ffm in their recent extension of the FEEC error estimates to semi-linear elliptic 
problems. Hence, it serves as an obvious assumption for extending our evolution results 
to the semi-linear case. 

For c M", define the Bochner semi-linear mixed weak form parabolic: Given / 
and g, find (m, a) e HT' x HT'-^ such that 

(Mi,0)-(diva,0) + (F(M),0) =(/,0), V</>Gi/A", tel, 

(cT, w) + (n, div =0, yujeHA''-\ tel, (7.3) 

uiO) =g. 

The semi-linear semi-discrete Bochner parabolic problem is thus: Find {uh, ah) e 

X A;;-^ such that 

{uh,u(j)h)-{divah,(ph) + {F{uh),(t>h) ={f,(ph), V0,,gA^, tel, 

{ah,ujh) + {uh,divuh) =0, Vw/,gA;J-\ tel, (7.4) 

where gh e A^ is an approximation of g. Analogously to the linear case, for any ^ L 
define the time-ignorant linear discrete elliptic problem: find {uh, (Jh) e k^x A^^ 
such that 

(diva;,,^;,) - (Az/(to),0h) =0, V 0/, G A^J, tel, 

+ (w/^divw^ =0, ^ ujh e kl'^ , tel, (7.5) 
w/i(0) = gh, 

where now u is the solution to the continuous semi-linear problem (|7.3I) . Similarly, define 

p{t) := Uh{t) - u{t), 

e{t) ■.= Uh{t)-Uh{t), 

e{t) := ah{t) - ah{t), 

where u, a and their discrete counterparts are now solutions to the corresponding semi- 
linear problems. We have an analogous lemma. 

Lemma 7.1. The semi-linear error functions satisfy the semi-discrete formulation, i.e. 

i9t,(j)h)- (dive, (j)h) =-{ptAh), ^(pheAl, tel, 
(e, Uh) + {0, div Uh) =0, ^ Uh e A^'^ tel. 
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Proof. The second equation is immediate from the second equations in (17.41) and (I7.5I ). 
The first equation can be written out as 

{uh,u <Ph) - (div ah, 4>h) + (div at, (ph) - {uh,t, (j^h) = {ut, 4>h) - {uh,t, M 
which is reduced as follows: 

iuh,u <Ph) - (div Oh, (ph) + (div dh, (ph) = (ut, <ph) cancel like terms 

iuh,t, (ph) - (div ah, (ph) = -(Am, (ph) + (ut, (ph) by (|73]) 

(/, (Ph) - {F{uh),(ph) = -{Au, (Ph) + {ut, (ph) by dH 

This says that the continuous problem Ut — Au + F{uh) = f should hold in a weak sense 
when tested against any of the functions in A^. This is guaranteed to be true since we 
chose C A" = L^. Thus, the error equations hold as stated. □ 

Theorem 7.2. Fix I := [0,T]. Suppose {u,a) is the solution to such that the 
regularity estimate 



holds for < s < Smax <^nd t E L Assume that the operator F satisfies the Lipschitz 
assumption ([7.2\) . Choose finite element spaces 



A n-l 
^^h 



Vr+iA^-\r) 



or 

1 n— 1 / 



n+iA"-i(r) 



= Pr+iA"(r) (= VrA^T)) 



Then for < s < Sma.x, dh defined by ([5.4\l , and {uh, ah) the solution to ([7.4\l , the 
following error estimates hold: 



ch 



\llh til |/^2^n ^ 



L2^i^L2) -r yT\\ft\\mi,L2)) ifr = 

ch'^' {\\f\yii,H^) + ^\\ft\\LHi,Hn) /'^'"^ >0,ifs<r-l 

(7.8) 



\(^h ~ (^\\£^2^n-l ^ < 



i \\f\\L^(l,H^) + V^II/t|lL2{/,L2)) 

ifr = 0,s = 0, Al-^ = A"-i(r) 
ifr = 0, s < 1, A^-i = PiA"-i(r) 
for r>0, ifs<r — 1 



(7.9) 



\di\{ah - (t) 



L2(I,H^) + h ||/t|lL2(/,L2; 



< < 



ifr = 0,s<l 

for r > 0, if s < r — 1 



(7.10) 
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Proof. The proof is very similar to that of Theorem 15 .21 Equation (17.51 ) is the k = n case 
of the discrete mixed variational problem examined by Hoist and Stern in [fTTl Equation 
(9)]. Therefore, we can use the same type of triangle inequality from (15.111) to recover 
the estimates. By [|17> Theorem 4.2], we have the estimates 

I ch''^^ \\f{t)\\^s if s < r - 1, for r > 




n+iA"-i(r) 



Mt)-a,{mL^<ch'^'\\m\\H^ 

(7.12) 

||div(a(t) - < ch'^Wfm^, , if s < r + 1. (7.13) 

An explanation of how these estimates are derived from the results of fTT| is given in Ap- 
pendix |Al Note that these estimates are exactly the same as the corresponding estimates 
(15.121 ). (15.141 ) and (15.161 ) from the linear case. The proof then proceeds exactly as before 
since the rest of the argument does not appeal to the linearity of the problem at all. □ 

This approach seems likely to extend to semi-linear hyperbolic problems as well. 
Since the well-posedness of such problems is a significant issue in its own right, however, 
we do not consider such an approach in the present work. 

8. Concluding Remarks 

In this article, we have extended the Finite Element Exterior Calculus of Arnold, Falk, 
and Winther |[3l [2l for linear mixed variational problems to linear and semi-linear para- 
bolic and hyperbolic evolution systems. Both the parabolic and hyperbolic cases make 
strong use of the smoothed projection operators vr^, which are one of the most elabo- 
rate and delicate constructions in the FEEC framework. In the parabolic case, the use of 
the TT^ operators was hidden somewhat by the use of elliptic projection error estimates, 
proofs of which rely on properties of these operators. In the hyperbolic case, the proof 
techniques use these properties more explicitly. In any case, the formal treatment and 
generalization of these operators by Arnold, Falk and Winther can now be seen as a 
useful tool for the analysis of evolution problems as well as elliptic PDF. 

We have also seen in this article how the recent generalizations of the FEEC by Hoist 
and Stern |fT6l [TTl for semi-linear elliptic PDF can be extended to evolution PDF as 
well, both parabolic and hyperbolic types. We also anticipate that the basic approach to 
analyzing variational crimes in |fT6l [TTl for the linear and semilinar elliptic cases will 
also work in the case of evolution problems; we will explore the question of variational 
crimes in a subsequent article, with the target being the analysis of surface finite element 
methods for evolution problems. 
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Appendix A. Explanation of Semi-Linear Error Estimates 

In this appendix, we explain why estimates (17.1 II ). (17.121 ). and (17.131 ) follow from [fTTl 
Theorem 4.2]. We will focus just on the r > case of (17.1 II) as it requires the sharpening 
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of a special case of an estimate appearing in [fTTl Theorem 4.2]. The other cases work 
out along similar lines by a direct application of the Hoist and Stern estimates. 

First, we recall some notation from [3] used in (VJ). If (W, d) is a Hilbert complex with 
associated domain complex (V, d) and parametrized subcomplex family (V^, d), denote 
the best approximation in PF-norm by 

E{w) = inf - , weW^. 

The relevant result from Hoist and Stern [fTTl Theorem 4.2] is stated as 

\\u - UhWy + \ \p-Ph\\w < c{E{u) + E{du) + E{p) 

+ r][E{a) + E{da)] + {5 + ^)E{da) + ^E{P^u)), 

(A.l) 

where rf, 6, and /i are coefficients defined as the norms of certain abstract operators, 
u G Wk, and p is a harmonic A;-form with discrete counterpart ph introduced to make the 
abstract Hodge-Laplacian problem well-posed. 

Casting this into the context of the deRham complex, we have 

{W, d) = (L^A, d) and (V, d) = {HA, d). 

Since we are interested here only in the case k = n, there are no harmonic A;-forms so 
that p = Ph = 0. Further, du = since rfA" = 0, whereby | |m — -u/il |y = | |m — -u/il = 
I |m — tt/i| 1^2- This eliminates the error terms in p and du, giving us the reduced estimate 

\\u - Uh\\L2 < c{E{u) + rj[E{a) + E{da)] + {5 + fi)E{da) + fiE{P^u)). 

Crucially, this estimate can be reduced further when k = n. The derivation of (|A.1|) uses 
the estimate 

\\d{u - Uh)\\w < c{E{du) + 7][E{da) + E{p)] 
from [[3l Theorem 2.11] which is unnecessary here since the left side is always zero. 
Since this is the only part of the derivation that requires the term r]E(da), we can drop 
it, yielding 

\\u-Uh\\L2 < c{E{u) +7]E{o-) + {6 + fi)E{d(r) + ^,E{P^u)). {A.l) 

We now give bounds on each of the terms in (IA.2I) . The coefficients appearing in the 
abstract estimates can be stated in terms of powers of h in the deRham context. These 
appear in [3^ p. 312] as 

= 0{h), 6 = 0(/i"^'"(2,r+i))^ ^ ^ 0{K'+^). 

To bound the error terms, Arnold, Falk and Winther define smooth projection operators 
TT^ : Li^A^ip) — 7- satisfying optimal convergence rates as stated precisely in [[3l 
Theorem 5.9]. For instance, if A\ is one of V^+iA'^iTh) or, if r > 1, VrA^iTh) then 

Ik-'^^^IL^AfeCf^) ^ c/i'll^ll//»AHn)' ioxw eH'A^{Q), 0<s<r + l. 

These types of results bound E{w) in terms of | Itf^l Iji^^^^fc, which is in turn bounded in 
terms of 1 1/| Ij^^^^Afe by the regularity hypothesis (17.71) . Summarizing these results, we have 

E{u)<ch^^'\\f\U. 

E{a)<ch^^^\\f%. 

E{da)<ch^\\f\\^. 

E[P^u)<ch^^^\\f%. 
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We can now prove (17.1 II) by collecting results and applying them to (IA.2I) . yielding 

The greatest common factor from the above expression is h^^"^ hence this is the overall 
order estimate that can be inferred, as was claimed. 
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